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Oh Abstract : 

Qh We calculate the functions of parton fragmentation into 3 Pj quarkonium at order a 2 , 

4^ where the parton can be a heavy or light quark. The obtained functions explicitly satisfy 

the Altarelli-Parisi equation and they are divergent, behaving as z~ x near z = 0. However, 
if one choses the renormalization scale as twice of the heavy quark mass, the fragmentation 
functions are regular over the whole range of z. 



1. Introduction 



Fragmentation functions are in general nonperturbative objects in the context of the 
factorization theorem in QCD[1]. This makes it hard to study them by starting directly 
from QCD. However, if partons, i.e., quarks and glouns, undergo fragmentation into a 
quarkonium, fragmentation functions can be factorized — they are sums of products of 
constants and coefficient functions. The constants represent the nonperturbative effects 
and may be defined as matrix elements in nonrelativestic QCD(NRQCD) or are related to 
a wavefunction, while the coefficients functions can be calculated with perturbative theory. 
With this factorization, fragmenation functions for various quarkonia were studied in [2-9]. 

A quarkonium contains a heavy quark Q and its antiquark Q } which move with a 
small velocity v in the quarkonium rest frame. Because of the small velocity the bound 
state effect, i.e., the nonperturbative effect in the quarkonium, can be well described by 
employing NRQCD. Recently, such an approach has been established[10]. The approach 
is basically distinct than early treatment within the color-singlet model(for a nice review 
see [11]). In the color-singlet model one treats a quarkonium system simply as a bound 
state of Q and Q } where the QQ pair is in a color-singlet state, and the nonperturbative 
effect is contained in the wavefunction of the bound state. Then an expansion in the 
small parameter v is made and to leading order only the wavefunction at the origin or 
its derivative at the origin is involved. This model has serious problems. First, it does 
not tell us how to handle the Coulomb singularities. Since an expansion in v is made, 
this type of singularity must appear because of massless photons and massless gluons. 
The effect related to the Coulomb singularities is nonperturbative. In the color-singlet 
model these singularities were absorbed into the wavefunction without solide reasons from 
theories. Second, infrared(I.R.) singularities appear when a P-wave quarkonium is involved 
and appear even at the leading order of a s . In the model the I.R. singularities were 
regularized as divergences in the limit of the zero-binding energy. Such I.R. singularities 
clearly indicate that the wavefunction of a P-wave quarkonium can not contain all the 
nonperturbative effects. In the approach [10], quarkonium systems are analysed in the 
framework of QCD. A systematice expansion in v can be made and the nonperturbative 
effect is represented by matrix elements in NRQCD. With these matrix elements one can 
show that Coloumb singularities are factorized into these matrix elements. From the point 
of view of a relativistice quantum theory a quarkonium system consists not only of a QQ 
color-singlet but also many other components like \QQG > ■ ■ ■ etc.. In the approach 
of [10] the effect of these other components is also included in the systmatic expansion 
in v. It should be emphasized [10,12] that a QQ color-octet state is as important as a 
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QQ color-singlet for a P-wave quarkonium. Therefore one should take both states into 
account. It is already shown that results at the one-loop level for gluon fragmenation into 
3 Pj quarkonium[9] are free from I.R. singularities and from Coloumb singularities. 

In this work we will study the quark fragmentation into a 3 Pj QQ quarkonium at 
leading order a s . We will show that because of the contribution of the color-octet QQ 
pair, a light quark q can also fragment into the quarkonium at the same order of a s as a 
heavy quark Q. The heavy quark fragmentation into P-wave quarkonium was originally 
studied in[8]. As pointed out in [9] the results for charmonium can not be correct, because 
the fragmenation functions obtained there do not satisfy the Altarelli-Parisi equation. 

The work is organized as follows: in Sect. 2 we introduce the defmition[13] of renormal- 
ized quark fragmentation functions and the factorzation[10] forms for quark fragmentation 
into 3 Pj quarkonium. Further details may be found in [10,13]. In Sect. 3 we start from 
the definiton to calculate the heavy quark fragmentation functions. In Sect. 4 we calculate 
the light quark fragmentation function. Finally we discuss and summarise the results of 
our work in Sect. 5. Throughout this work we always assume that the polarization of the 
quarkonium is not observed. We will use dimensional regularization and work at leading 
order in the expansion of v. 



2. The Definition of Quark Fragmentation Function 
and the Factorization Form for 3 Pj quarkonium 



As we will use dimensional regularization we give the definiton of the quark fragmen- 
tation function in d dimensions. To give the definitions for a fragmentaion function it is 
convenient to work in the light-cone coordinate system. In this coordinate system a (/-vector 
p is expressed as p^ = Pt), with p + = (p° + p d ~ 1 ) / \/2, p~ = (p° — p d ~ 1 ) / \/2. 

Introducing a vector n with n' 4 = (0, 1, 0, • • • , 0) = (0,1, Ot), the fragmentation function 
for a spinless hadron H or for a hadron without observing its polarization is defined as[13]: 

D( H/ q ( z )=^- J ^"e" jP+x " / ^Tr color ^Tr Dirac {n- 7 <0| (Z (0) 

r-OO 

Pexp{-ig s / d\n ■ G T (An' t )}a{ f (P + , T )a ff (P+, T ) (2.1) 
Jo 

poo 

Pexp{ig s / d\n ■ G T (An' t )}g(0, , O T )|0 >, 

J x~ 

where G^x) = G^(x)T a , G^(x) is the gluon field and T a (a = 1, • • • , 8) are the SU (3)-coIot 
matrices. The subscript T denotes the transpose. aJj(P) is the creation operator for the 
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hadron H. For hadrons with nonzero spin the summation over the spin is understood. 
The definition is a unrenormalized version. Ultraviolet divergences will appear in D^ q (z) 
and call for renormaization. Following [13] the renormalized gluon fragmenation function 
can be defined as: 



D H/Q (z) = D%) Q {z) + f J L ^ y -) D H/a(y)- (2-2) 

a = G,q z 



Here the summation over all possible partons is understood. The function L a (z) is chosen 
so as to cancel the U.V. divergences. In MS scheme L a (z) takes the form: 

L a (z) = J2^L^ } (zl (2.3) 

N 

where e = 4 — d. From Eq.(2.2) one can derive the Altarelli-Parisi type evolution equation 
for the fragmentation function. We will use the modified MS scheme, where L a (z) is chosen 
to cancel the terms with N e = -| — 7 + ln(47r) . The function D}j/ q (z) is interpreted as the 
probability of a quark q with momentum k to decay into the hadron H with momentum 
component P + = zk + , it is gauge invariant by definiton. Further, it is also invariant under 
a Lorentz boost along the moving direction of the hadron and under a rotation with the 
direction as the rotate axis. 

If the hadron is a 3 Pj quarkonium, a factorized form for the fragmentation function 
can be taken. We will use the notation \J f° r the 3 Pj quarkonium. At the leading order 
of v Dnj q {z) can be written according to [10] as: 

D xj/q (z) = < 010*^)10 > +^ < 010*^)10 > . (2.4) 

where Di(z^J) and D$(z) are dimensionless and D$(z) is same for all J. The operators 
Of J ( 3 Pj) and Og J ( 3 Pi) are given by: 

Of ( 3 Si) = x^ t T a 4,(a^ H a H )^a t T a x 

Of ( 3 Po) = IxH, ~ \ 6 .a)^ H a H )^{ - 1 6 .*) x 

Of ( 3 Pi) = \x\ - \ 6 xa)^ H a H )^( - \ 6 x,),X ^ 
Of ( 3 P 2 ) = xH ~ \ D {1 a j} )^(a^a H )^( - % - D {l a j} ) x , 

where D is the space part of the covariant derivative D^ 1 and ai(i = 1,2,3) is the Pauli 
matrix. The notation {ij} means only the symmetric and traceless part of a tensor taken. 
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In Eq. (2.5) tp and are fields with two components for the heavy quark Q and its 
antiquark Q in NRQCD. M is the mass of the heavy quark. d\j is the creation operator 
for the hadron in its rest frame. The matrix elements in Eq.(2.4) are defined in NRQCD. 
In Eq.(2.4) the part with Dg is the contribution from a color-octet QQ pair in a 3 Si 
state and the part with D\(z, J) is the contribution from a color-singlet QQ pair in a 3 Pj 
state. We will call them as the color-octet and color-singlet components respectively. The 
matrix elements represent the nonperturbative effect, while D\ and Dg can be calculated 
perturbatively and they should be free from I.R. singularities. 

A good method to calculate D\ and Dg is to use wavefunctions to project out different 
states from a general QQ pair. At the leading order of v the projection can easily be worked 
out, details may be found in [14]. We will use a radial wavefunction Ri{r) to project the 
3 Pj color-singlet QQ state and an octet radial wavefunction R% a \r) to project the 3 Si 
color-octet QQ state. Calculating with these wavefunctions the l.h.s of Eq.(2.4) and the 
matrix elements in the r.h.s of Eq.(2.4) we can extract the functions D\(z, J) and Dg(z). 
The results for D\(z, J) and Dg(z) are independent on these wavefunctions. At the order 
of a s we consider, only the tree-level results for the matrix elements are needed, they are: 

< o\o*>CPj)\o > = 9(2 ^ + 1 V i(o)l 2 , 
<o|o 8 ^( 3 5 1 )|o> = f £|4 C) (0)I 2 , 

07T z — ' 

c 

where R'i(0) is the first derivative of i?i(r) at the origin. 

3. The Heavy Quark Fragmentation Function 

From the defintion in Eq.(2.1) we can always decompose the fragmentation function 
by sandwiching the operator \X >< X\ between ajj and ajj as: 

</q(-)~5>{--74M, (3-i) 

X 

where Th may be called the fragmentation amplitude for Q — > H + X . Here we have the 
conservation of total momentum only in the +- direction. 

3.1 The Color-Singlet Component 

The color-singlet component receives nonzero contributions at order a 2 s . The Feynman 
diagrams for Th are given in Fig.l. Because the QQ pair is in a color-singlet state, 
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there are two gluon lines attached to the quark line. Here, there are no divergences, so 
renormalization is not required. The calculation is complicated. Because of the summation 
over intermediate states we encountered integrals of the type: 



( dq T d -2 ( 2 (2-zfM 2 n no AK (o \ 

) y q T + 2 ) , for n = 2,3,4,5. (3.2) 

Z7T Z Z 



Here qx is the transversel momentum of the quark as the intermediate state in Fig.l. The 
integrals are finite and after performing the integrations we can extract: 

a / ~t , 16 , z(l — z) 2 



729 syrj (2-zf 
(192 + 384z + 528/ - 1376/ + 1060/ - -376/ + 59/), 



a / t x 64 ,. . z{\ — z) 
D 1 (z,J = l) = —a 2 3 ( t i) K ) 



729 sKrj (2-zf (3.3) 
• (96 - 288^ + 496/ - 408/ + 202/ - 54/ + 7/), 

A , T x 128 o, s z(l - z) 2 

DUz, J = 2) = a Ju,)-? ^ 

V ' ; 3645 sKHJ (2-zf 

■ (48 - 192z + 480z 2 - 668z 3 + 541/ - 184/ + 23/). 

These results agree with those in [8]. It is interesting to note that there is a common factor 
z(l — z) 2 (2 — z)~ 8 for all J, whereas there is a common factor z(l — z) 2 (2 — z)~ 6 for heavy 
quark fragmentation into S"-wave quarkonium[6]. Note that the same diagrams in Fig.l 
contribute to heavy quark fragmentation into S"-wave quarkonium. The difference between 
these two factors are because for P-wave quarkonium the derivative of the fragmentation 
amplitude with the relative momentum between Q and Q is involved whereas for S"-wave 
quarkonium only the amplitude itself is involved. The appearance of these common factors 
can be roughly understood by counting the denominators due to the quark- and gluon 
propagators in the amplitudes and factors from phase-space. A successful model for heavy 
quark fragmentation was obtained through this way [15]. However, we will see in the next 
subsection that such a counting rule will be violated due to renormalization. 

3.2 The Color-Octet Component 

For the color-octet componenet, not only the diagrams in Fig.l, but also the diagrams 
in Fig. 2 will contribute The two diagrams in Fig. 2 were missing in [8] and they lead to 
divergences. Instead of integrals in Eq.(3.2), we have: 

( dq T sd-2 ( 2 , (2-z)M n ( 
(-^jr) (<1t + j , for n = 1,2,3,4. (3.4) 
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The integral with re = 1 is ultraviolet divergent, requiring renormalization. For the renor- 
malization we note that the function of gluon fragmentation into \ j quarkonium is nonzero 
at order a s [9,6,7]: 

D xj/G (z) = ^a s S(l ~z)~< OlOrCS^O > +0(a 2 s ). (3.5) 

Substituting this into Eq.(2.2) for H = xj we can easily chose the function Lo(y) to cancel 
the divergence. Finally, we obtain the renormalized function Dg(z): 

M*) = ^M{^(1 + (1 - ,) 2 )(ln (J^) - 21n(l - Z -))-z 

2(1 -z) (3 - 6) 
+ -f ^ 192 ~ UMz + 2016 * 2 ~ 1360* 3 + 352z 4 - Uz 5 - 5z 6 )}. 

Here there is is no common factor like z(l - zf(l - z)~ 6 . With the results in Eq.(3.6) 
and in Eq.(3.3) we complete the heavy quark fragmentation function at order a 2 s . This 
function should in general satisfy its evolution equation: 



^ h/q( ,f*) = ^ I ^Lp Q ^ q ( z /y }fi )D H/q (y } fi)+ I —P Q ^ G (z/y,fi)D H / G (y,fi). 

Jz y Jz y 

(3.7) 

The spliting functions PQ^ g (y,/i) and Pq^g(j/, /") are in the one-loop approximation the 
same as those for parton distributions. Using this fact and the result in Eq.(3.5) we obtain 
the evolution equation for the quark fragmenation at order of a 2 s : 

= ±a%M± < 0|O 8 -(^)|0 > .i + li^l. (3.8, 

Substituting our results into the l.h.s. in Eq.(3.8) one can check that our results are in 
agreement with this equation. From our result the fragmentation function is divergent as 
z~ x when z — > 0. However, this singularity disappears if we chose the renormalization 
scale /i as twice of the mass M. The same was also found in the gluon fragmentation in [9]. 
This property is important for possible applications of our result. In practical applications 
one solves the evolution equation at /i numerically, where one needs the moments of our 
result for the fragmentation function at some initial scale /io as the boundary condition. 
To insure that the perturbative result is a good approximation, one should chose /io ~ M 
to avoid large logarithmc terms in higher order. Our result tells that one should chose 
/io = 2M to avoid these terms in higher order and to safely calculate the moments in other 
hand. 
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4. The Light Quark Fragmentation Function 



Since a color-octet QQ pair will lead to a contribution to P-wave quarkonium pro- 
duction at the leading order of v, a light quark q can undergo fragmentation into XJ by 
generating a color-octet QQ through emission of a virtual gluon. Such a process happens 
at the same order of a s as the heavy quark fragmentation. The Feynman diagrams for Tjj 
are those in Fig. 2, where the quark line attached by the double line is for the light quark 
q. At the leading order of a s and v the light quark fragmentation function D XJ / q (z) can 
be written: 

D xj/q (z) = ^fi- < 010^)10 > . (4.1) 

The color-singlet component only becomes nonzero in higher order than a 2 . The calcula- 
tion is similar to the previous section. We introduce the notation y = -^j-, where m q is the 
mass of q. The result for D$ jq (z) is: 



Ds, q (z) = + (1 " -) 2 )[ln (^) - In ((1 - z) + \y 2 z 2 ] 



z -2(l-z)(2 + y 2 ) 



2- z 



(4.2) 



4(1 - z) + y 2 z 2 



Again the light quark frgamentation function must satify its evolution equation. At order 
a 2 s this equation is the same in Eq.(3.8). It is easy to check that the function in Eq.(4.1) and 
(4.2) satisfies the evolution equation. The light quark fragmenation function has the same 
property near z = as the heavy one, i.e., it is divergent as z~ x at any remormalizaion 
scale /i except when /i = 2M. The light quark mass m q can be safely neglected. With 
m q = the function in (4.2) becomes: 

Ds, q (z) = ^o?M{\{\ + (1 - ,) 2 )[ln (J^) - ln(l - z)] - 2z). (4.3) 

For the convenience of later discussions we introduce here some relations between 
the various matrix elements in Eq.(2.5). In principle these matrix elements have series 
expansions in v and the leading order is v 2 . Since we only work at the leading order, the 
higher order corrections can be neglected. In this case, the matrix elements in Eq.(2.5) are 
related to each other with a spin factor of \j. We introduce two parameters H\ and H' s 
as in [10], and the relations can be expressed as: 

< 0\O^ j CPj)\0 >^(2J + 1)M 4 H u 

(4.4) 

< 0|O^ J ( 3 5i)|0 > w (2 J + 1)M 2 H' 8 . 
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With these relations the whole set of the quark fragmention functions contains only two 
unknown parameters, which can only be computed nonperturbatively or extracted from 
experiments. 



5. Discussion and Summary 



Some useful quantities of parton fragmentation functions are their first moments. 
These moments allow one to roughly estimate a single hadron production rate though 
fragmentation, where the rate may be taken as product of a parton production rate and 
the corresponding first moment, where the summation over different partons is understood. 
We will give results of the first moments of our fragmentation functions. We denote the 
first moment as M(qf — > xj)i where qf stands for Q or q. Taking /i = 2M, we obtain: 

M (q ^ XJ ) « 0.029(2J + hnVM/!^. 

M(Q ^ X o) ~ 0.024a*(2M)^ + 0.035^(2M)|l 

H' H (5-1) 

M(Q - xi) ~ 0.072a*(2M)-^ + 0.039^(2M)^, 

M(Q X2) ~ 0.12^(2M)^- + 0.015^(2M)^. 

Here we neglect the mass of the light quark. We take the c quark as an example to 
give some value for the moments. For the value of H\ and H' s we use the estimates in 
[12]. Taking M = m c = 1.5GeV and a s (2m c ) = 0.26, the value for M(q — > Xcj) is 
3.8 x 10~ 6 (2J + 1) and the value for M(Q -> X cj) is 2.6 x 10" 5 , 3.4 x 10" 5 and 2.5 x 10" 5 
for XcO, Xci and Xc2 repectively. The contribution from the color-octet componennt is not 
negligible. For Xc2 the contribution from the color-octet component is roughly 70% of 
the heavy quark moment. From these values one can see that the moments for the light 
quark fragmentation are roughly one order of magnitude smaller than those for the heavy 
quark. But the contribution from light quarks should not be neglected, especially in a 
hadron reaction, since the production rate of light quarks as parton may be larger than 
the production rate of a heavy quark and hence a substantial contribution from light quark 
fragmentation to the XJ production is possible. 

With the results here and those in [9] the functions of all possible parton fragmenation 
into 3 Pj quarkonium are calculated at order a 2 s . Only two parameters, which represent the 
nonperturbative effect at the leading order of v, are not known precisely. The functions 
have the general feature that they are divergent as z~ x when z — > 0. But at /i = 2M 



8 



they are regular distributions over the whole range of z. The functions also satify the 
Altarelli-Parisi equation, as expected. 

Acknowledgment: The author would like to thank Dr. A. Rawlinson for reading the 
text carefully. This work is supported by Australain Research Council. 
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Figure Captions 

Fig.l: The Feynman diagrams for the color-singlet component of the heavy quark frag- 
mentation. The line is for the heavy quark, the wavy line is for gluons. The double line 
represents the line operator in Eq.(2.1). 

Fig. 2 The Feynman diagrams for the color-octet component of the heavy and light quark 
fragmentation. 
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Fig. 2 



